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Abstract 

This paper presents techniques and analysis of time-of- 
arrival (TOA) based ranging with analog chip-radios used 
for IEEE 802.11 (Wi-Fi). TOA and round-trip-time estimation, 
reconstruction of a wideband signal, and local oscillator 
clock frequency-offset estimation are addressed. Theoretical 
performance analysis is performed and concise formulas are 
derived. Both the theoretical and experimental results dem- 
onstrate that by means of these techniques it is feasible to 
achieve sub-meter ranging accuracy in severe non-line-of- 
sight (NLOS) indoor environments using the same analog 
chip radios as used for Wi-Fi data communications. 
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Introduction 

The development of mobile radio systems has been 
dramatic, partly due to the development of chip radio 
systems, including both the analog radio transmitter 
and receiver. Initially, hardware was based on both 
analog radio chips and a separate digital signal proc- 
essing chip, but later versions integrated both these 
functions into a single chip. While the main focus of 
these developments has been on digital radio commu- 
nications, in more recent times the availability of inte- 
grated single chip Global Positioning System (GPS) 
receivers (Samper et al 2009, Diggelen 2009) has al- 
lowed the rapid development of position-based ser- 
vices in devices such as "smart" phones. Simultane- 
ously, radio technology has been extended to indoor 
environments for wireless data communications, par- 
ticularly the widespread adoption of IEEE 802.11 (Wi- 
Fi) as the defacto standard. With the availability of 
integrated chips, Wi-Fi has become standard in both 
Notebook/Tablet computers and in smart phones, with 
Wi-Fi base stations for data communications currently 
common in both workplace and home. Investigations 


on Wi-Fi positioning started more than one decade 
ago and since then Wi-Fi positioning has drawn sig- 
nificant attention mainly due to its simplicity. For in- 
stance, it has been used in commercial applications 
such as tracking mobile phones indoors when the GPS 
fails to work. However, accurate (say better than 
1 meter) radiolocation is not commonly available, al- 
though applications such as tracking people (security, 
first responder/emergency and sporting applications) 
and assets (Bahl & Padmanabhan 200, Fontana et al 
2003, Fontana 2004, Hedley et al 2008, Mourad et al 
2011) are commonly cited as future developments with 
such accuracy as a requirement. 

In indoor environments GPS may not work well since 
the satellite signals are blocked by the buildings, and 
typically non-line-of-sight (NFOS) multipath condi- 
tions predominate. Accuracy requirements for dense 
urban environments and indoor locations are typically 
much greater than the corresponding requirements for 
outdoor positioning. For example GPS accuracy out- 
doors is typically around 10 meters (Diggelen 2009), 
while for the indoor case with applications such as 
tracking people, an accuracy of 1 meter or better is de- 
sirable to allow location determination, such as to a 
particular room inside a building. Although GPS re- 
peaters may be used to obtain better performance 
(Jardak & Samama 2009), in general GPS cannot meet 
these operational requirements indoors, and thus 
some other technology is required. This paper investi- 
gates the possibility that the chip radio used for Wi-Fi 
data communications could also be used for indoor 
position determination with 1-meter accuracy. 

Because of the complexity of the indoor radio propa- 
gation environment, it is difficult to achieve 1 -meter 
radio ranging accuracy. While UWB systems can have 
good indoor ranging error performance of the order of 
10-20 centimeters, the severe limitation on the avail- 
able transmitter power restricts the range to typically 
10 meters (Roy et al 2004, Ghassemzadehl et al 2005, 
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Alavi & Pahlavan 2006, Bellusci et al 2008, Gentile & 
Kik 2006, Alsindi et al 2009). Thus, to cover large in- 
door spaces, a large number of nodes is required, 
which is logistically difficult and expensive. An alter- 
native is to use the industrial, scientific and medical 
(ISM) bands (particularly the 2.4 GHz and 5.8 GHz 
bands) as used by Wi-Fi devices, where the much 
higher allowable transmitter power (FCC Part 15) re- 
sults in indoor ranges of typically 40-80 meters, and 
hence far fewer nodes are required to cover the same 
area. However, with maximum bandwidths of respec- 
tively 80 MHz and 150 MHz, the time resolution of 
ISM-band positioning systems is much less than that 
of UWB systems, and the positional accuracy is con- 
siderable worse (Yu et al 2009). While the development 
of specific new technologies may meet future needs, a 
commercially attractive alternative is to adapt the cur- 
rent Wi-Fi data communications technology for high- 
accuracy radiolocation. It is noted that the 20 MHz Wi- 
Fi channel bandwidth will have a nominal time resolu- 
tion (pulse rise-time) of the order of 100 ns (or about 
30 m), but the actual ranging accuracy depends on 
many other factors including the operating environ- 
ment (outdoors/indoors), signal-to-noise ratio, quality 
of the analog radio hardware, and the signal process- 
ing algorithms. Section 2 provides an overview on ac- 
tual performance of some Wi-Fi based ranging sys- 
tems, which typically have an accuracy of 2.5 m to 5 m. 
As the Wi-Fi radio has not been designed for this pur- 
pose, various aspects and limitations must be over- 
come for accurate TOA measurements to the order of 
1 ns. Such techniques are described in detail in the fol- 
lowing sections, including performance estimates, and 
confirmation by actual measurements using a chip 
radio (Maxim MAX2828). While this equipment is ex- 
perimental in nature, the results show that future inte- 
gration into Wi-Fi chips could provide the dual roll of 
data communications and position location. It is also 
worth noting that one of the recent research focuses in 
Wi-Fi positioning is on the integration of such a tech- 
nology with other technologies such as GNSS, RFID, 
or Dead-reckoning (Retscher & Fu 2009, Li et al 2010). 

The main contributions of this paper can be summa- 
rized as follows: 

1) Propose a technique for concatenating multiple 
Wi-Fi frequency channel data to effectively ob- 
tain a wideband signal, and derive theoretical 
formulas to estimate the degrading in perform- 
ance associated with this technique compared 
with the wideband signal case. 


2) Develop a procedure for the correction of local 
oscillator frequency offsets, and show that an 
accuracy of the frequency corrections of the or- 
der of 1 ppb can be achieved for mobile nodes. 
Previous methods described are only valid for 
the stationary case (anchor nodes). 

The remainder of the paper is organized as follows. 
Section 2 provides an overview on related work. Sec- 
tion 3 reviews a technique of range determination us- 
ing time-of-arrival (TOA) data without the need for 
time synchronization, with particular application to 
Wi-Fi radios. Section 4 describes and analyzes the per- 
formance of a method of using Wi-Fi 20 MHz channels 
to achieve an effective wideband channel, and hence 
greater time resolution for TOA measurements. Sec- 
tion 5 describes and analyzes a technique for very ac- 
curate frequency synchronization throughout a mesh 
network in an indoor NLOS propagation environment 
with stationary and mobile nodes. Section 6 provides 
some measured performance results using the meth- 
ods described in previous sections. Finally, Section 7 
summarizes and concludes the paper. 

Overview of Related Work 

Before considering the specific details of the proposed 
approach, this section briefly overviews some related 
work in the general area of Wi-Fi positioning. 

1) Intel research has developed a Wi-Fi position- 
ing system based on TOA ranging, which re- 
quires some hardware modifications to be im- 
plemented (Golden & Bateman 2007). The 
range is computed through performing a num- 
ber of different exchanges of messages between 
an access point and a stationary mobile. Al- 
though some advanced signal processing tech- 
niques are used, the range root-mean-square 
error (RMSE) associated with 5 access points is 
reported as 2.66 m. Performance data for a 
moving node is not provided. 

2) Typically in currently proposed systems, two- 
way TOA (equivalent to a single round-trip- 
time (RTT)) measurements are employed for 
range measurements, such as using the packet 
pairs DATA and ACK in the IEEE 802.11 stan- 
dards. In (Llombart et al 2008) a purely soft- 
ware-based ranging and positioning approach 
is proposed and implemented by using the ex- 
isting Wi-Fi hardware. Instead of using packet 
pairs, packet sequences (RTS, CTS, DATA and 
ACK) are exploited to generate multiple RTT 
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measurements, resulting in improved meas- 
urement accuracy, although still only 4 m. 

3) The development of a WLAN channel sounder 
for IEEE 802.11b is reported in (Ciurana et al 
2009). The sounder directly uses the physical 
layer characteristics of Wi-Fi systems to esti- 
mate the channel impulse response. Using a 
multipath propagation model and conducting 
a Monte-Carlo simulation, the resolution of the 
system is found to be nine nanoseconds, 
equivalent to 2.7 m. 

4) A strategy is proposed in (Jemai & Kumer 2008) 
to reduce the position determination latency 
without degrading positioning accuracy. In 
particular a method is described to optimize 
the tradeoff between a short RTT response time 
and position update period. (Hoene & Will- 
mann 2008) investigates the scalability issue as- 
sociated with the Wi-Fi network based posi- 
tioning. 

None of the above is related to actual working posi- 
tioning systems, and thus the actual performance in a 
real operating environment cannot be assessed. To 
investigate the potential of using Wi-Fi chips for posi- 
tion determination to achieve high-accuracy position- 
ing, the techniques described in this paper have been 
applied to the design of an experimental positioning 
system called the Wireless Ad-hoc System for Posi- 
tioning (WASP) (Hedley et al 2008, Sathyan et al 2011, 
Sharp, Yu & Hedley 2012, Sharp, Yu & Sathyan 2012). 
The system using a commercially available Maxim 
MAX2828 Wi-Fi chip in each node in the network op- 
erates in the 5.8 GHz ISM band with an effective 
bandwidth of 125 MHz synthesized from 20 MHz 
channels, and has been designed to achieve indoor 
positioning with an accuracy of 1 m (or better) based 
on a mesh network of nodes. While the above- 
referenced WASP papers give some information on its 
design and performance, key aspects such as the gen- 
eration and detection of the wideband signal and the 
frequency synchronisation in the important case of 
mobile nodes have not been described previously. 
These topics are the main focus of this paper, provid- 
ing further information on how the performance of the 
WASP system is achieved. 

Range Measurement Overview 

This section reviews the general theory of determining 
ranges between nodes in a network using TOA meas- 
urements. Although some of the material in this sec- 


tion has been described elsewhere (Hedley et al 2008, 
Sathyan et al 2011, Yu, Guo & Hedley 2009), the fol- 
lowing provides more in-depth information related to 
the use of chip radios, and also provides the underly- 
ing theoretical framework for subsequent sections, 
particularly for the important practical case of a mo- 
bile node. 

Transmitter to Receiver Timing Measurements 

The determination of the propagation delay used for 
position estimation can be based on the measurement 
of the TOA of a radio signal at a receiver. However, 
the TOA measurement encompasses the total trans- 
mission delays from the transmitter to the receiver as 
measured relative to the clocks in the (mobile) trans- 
mitter and the receiver (stationary node or anchor 
node/base station), and includes the delays in the 
baseband hardware, the RF and IF sections of the ra- 
dios, as well as the propagation delay. A block dia- 
gram of the transmission path and the associated pa- 
rameters is summarized in Fig. 1. Note that transmis- 
sion can be from node m to node n or in the reverse 
direction. In this initial analysis, it will be assumed 
that the clocks in the transmitter and receiver have the 
same frequency, but later this restriction will be re- 
moved. However, it is important to note that the 
clocks are not synchronised in time, so that the clocks 
will have a random time offset relative to an absolute 
time reference. 



FIG. 1 BLOCK DIAGRAM SHOWING THE MAJOR COMPONENTS 
IN THE TRANSMISSION PATH FROM THE TRANSMITTER TO 
THE RECEIVER. 

The transmitter baseband electronics has a clock (local 
oscillator) which has an (unknown) time offset relative 
to absolute time. Further, it is assumed that the trans- 
missions are based on a pseudo-random noise code (or 
pn-code), which repeats at the code period; and that 
the code length is much greater than the total propaga- 
tion delay, so code ambiguity does not occur. As the 
time measurement has a repetition period, it is appro- 
priate for the (absolute) clock offset to be referred as a 
phase offset (</>) rather than an (unknown) absolute 
time offset. The baseband pn-code is used to modulate 
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the RF signal in the radio, which outputs the radio 
signal from an antenna. This process involves delays 
in the radio, filters, and cables, all of which are of un- 
known value. The total delay in the transmitter proc- 
ess is defined as A f * for node m. Note that the delays 
in the radio transmitter, particularly the filters can be 
several hundred nanoseconds - see Fig. 2, considerably 
greater than the propagation delay between the nodes, 
and further these hardware-related delays vary over 
time. The reason for the variability is mostly associ- 
ated with filter delay variations due to temperature 
and ageing effects. Further the length of the cable from 
the radio to the antenna can vary from node to node. 

Now consider the TOA measurement in more detail. 
The basic time period in the transmitter and the re- 
ceiver is the length of the pn-code symbol period (r ). 

As the code repeats with this period, normal continu- 
ous time is not available in the receiver, so that the 
transmitter and receiver time is related to absolute 
time by 

([> = t mod T pn 

If the clocks are assumed to be accurately synchro- 
nised in frequency (see Section 5 for details of the fre- 
quency synchronisation), then the local time in both 
the transmitter (tx) and the receiver (rx) can be defined 
by a phase offset determined at local clock time (/> = 0 . 
With reference to Fig. 1, the TOA measurement from 
node m to n (and ignoring measurement noise) can be 
expressed as 

TOA mn = (j> m + A“ + R mn lc + A" - </> n (2) 

where (/> is the clock offset phase (0 <tj>< T pn ) and c is 
the propagation speed of 0.2998 m per nanosecond. 



FIG. 2 MEASURED VARIATION IN THE DELAY PARAMETER 
FOR THREE DIFFERENT HARDWARE MODULES AS A 
FUNCTION OF TIME. THE DELAY DATA ARE MEASURED 
EVERY SECOND AND AVERAGED OVER A PERIOD OF THREE 
MINUTES FOR EACH RUN (RUN NUMBER) OF THE 
EXPERIMENT. THE PARTICULAR CHIP RADIO IS A MAXIM 
MAX2828. 


Thus the TOA measurement involves two parameters 
in the transmitter and two in the receiver, but as each 
node has both transmit and receive functions, each 
unit has three parameters associated with it (the clock 
phase, the transmitter delay, and the receiver delay). 
While (2) is perfectly satisfactory for determining the 
TOA and further data processing, some simplification 
is possible by the adoption of a different time refer- 
ence in each module. In particular, if the clock refer- 
ence is altered so that the timing reference point is at 
the transmitter antenna (rather than at the baseband 
clock), then the TOA expression becomes 

TOA mn = (, 1> m - A* ) + A* + R mn /c + A” - (*„ - A* ) 

= </ m +R mn / c + ( A n + Ki) ~ &n = Qm + R mn I c +A„ ~ 

( 3 ) 

Thus the TOA measurement can be expressed in terms 
of the clock phases in the transmitter and the receiver, 
the propagation delay, and the combined transmitter 
and receiver delays in the receiving unit only. As a con- 
sequence, only two parameters in each node need to 
be defined, namely the local clock phase and the com- 
bined radio delays. As the clock phases and radio de- 
lays are unknowns to be determined as part of the po- 
sition determination process, the processing is some- 
what simplified, as only two (nuisance) parameters 
need to be measured instead of three as in the original 
TOA expression. Therefore, (3) is the basic starting 
point in position determination. 

For typical positioning systems, nodes both transmit 
and receive, so that it is useful to consider the TOA 
measurements between a pair of nodes. Thus applying 
(3) twice, the resulting two TOA measurements are 

TOA mn =</>,„ +R mn /c + A n -<t>„ 

TOA nm = (j) n + R nm lc + A m - (/> m (4) 

RTT mn = ( T °A m „ + TOA nm )/ 2 = R mn /c + (A m + A n )/2 

where reciprocity in the propagation delay has been 
used (Jordan & Balmain 1968). Thus by adding two 
TOA measurements into a single RTT measurement, 
the clock phase parameters can be eliminated. If the 
radio delay parameters are known (preset) values, 
then the summed TOA measurement can be used to 
estimate the range. However, this typically results in 
low accuracy measurements due to the variability in 
the radio delays over time (Fig. 2), and will require 
individual calibration for greater accuracy. Such cali- 
bration techniques are beyond the scope of this paper, 
but described in (Hedley et al 2008, Sathyan et al 2011) 
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A key feature in the formulation described above is the 
use of pn-code phase rather than time. It is important 
to note that (4) is based on clock phase but not time, so 
that the two transmissions can be separated in time by 
a considerable amount, and thus an immediate ac- 
knowledgment of a transmission is not required. This 
subtle difference has important ramifications on the 
design of a mesh network. In a typical implementation 
of RTT, the receiving node replies immediately to the 
transmitting node. Therefore, in a network of N nodes 
a total on N(N - 1) time slots are required to avoid 
mutual interference; for even modestly sized mesh 
networks this arrangement is impractical as the update 
rate is low due to the large number of time slots. How- 
ever, as the RTT reply transmission can be delayed in 
time if the measurement is based on the pn-code 
phase as in (4), the N(N - 1)/2 ranges in the mesh 
network can be determined with only N time slots, in 
conjunction with the corresponding mesh network 
data communications. However, this method requires 
that the frequency offsets in the local clocks in each 
node is accurately compensated for, as (say) a typical 
local oscillator with 1-ppm accuracy would result in a 
RTT error of 1000 ns (300 meters) with a 1 second de- 
lay. Further, if the node moves significantly during the 
acknowledgment period, further timing errors will 
occur. 

Timing Measurement Errors 

The analysis in the previous section assumed ideal 
error-free TOA measurements, but in an actual system 
there will be measurement errors associated with mul- 
tipath propagation and receiver noise. To minimise 
multipath errors, TOA measurements need to be based 
on the leading edge of the despread received wide- 
band spread-spectrum signal. For the WASP system, a 
leading edge threshold algorithm (Sathyan et al 2011) 
is used, so the following discussion will be limited to 
that particular implementation; other leading edge 
algorithms (Yu, Sharp & Guo 2009, Sharp, Yu & Guo 
2009) have similar characteristics. 

The theoretical performance of the threshold algo- 
rithm in the presence of Gaussian noise can be calcu- 
lated by assuming that the nominal pulse shape (auto- 
correlation function of the pn-code) is triangular, cor- 
rupted with random noise with zero mean and a STD 

2 

such that the signal-to-noise ratio is y = l/cr„“, where 
the amplitude (A) of the pulse without noise is nor- 


malized to unity. The threshold for the detection of the 
leading edge has a nominal amplitude of a th A , or 
simply a lh in the normalized form. Although the tri- 
angular shape (with a rise-time of z pulse ) is only an 

approximation for the bandlimited signal case, pro- 
vided the threshold level is not too small (a lh >0.05) , 
typical shape distortions are sufficiently small (Sharp, 
Yu & Guo 2009) so that they have a minor effect on the 
TOA error due to Gaussian noise. Using techniques 
similar to those described in (Sharp, Yu & Guo 2009) it 
can be shown that the STD in the TOA error measured 
by the threshold algorithm is given by 


a 


£ 




( 5 ) 


where the last approximation is made because the 
threshold amplitude is small in order to reduce the 
effects of multipath signals. This result is used later in 
Section 5 in the analysis of the frequency synchroniza- 
tion performance. 


Multi-segment Concatenation 

To realize the design goals outlined in Section 1 the 
effective time resolution of the TOA measurements 
must be of the order of 10 ns, or alternatively a band- 
width of at least 100 MHz. While such a bandwidth is 
available in the 5.8 MHz ISM band, Wi-Fi chip radios 
designed for IEEE 802.11a operation have a bandwidth 
of 20 MHz, and thus it cannot meet this design re- 
quirement. The transmission approach adopted in the 
WASP system is for the transmitter to sequentially 
scan across the ISM band one channel at a time, and 
for the receiver to reconstruct the wideband signal by 
appropriate signal processing of the logged channel 
data. However, as the receiver must be retuned for 
each channel, the data from each channel are not glob- 
ally coherent (although coherent within each channel), 
even after an appropriate settling period for the re- 
ceiver local oscillator; a coherent wideband signal is 
essential to perform the despreading of the wideband 
spread-spectrum signal. One possible method to en- 
sure a coherent total signal is to overlap the bands in 
frequency, so that the overlapping frequency sections 
can be used to determine the differential phase be- 
tween adjacent channels. In the WASP system, eight 
overlapping bands of bandwidth 18 MHz (total 
144 MHz) are concatenated into a 125 MHz wideband 
signal, with about 2.5 MHz of overlap between adja- 
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cent channels. 

Procedure for Transmitting/Reconstructing the Wide- 
band Signal 

As this paper focuses on the complications associated 
with the use of a chip radio rather than the details of 
signal processing, the following is only a summary of 
the procedure to generate the multi-segment signals 
for transmission, and the subsequent reconstruction in 
the receiver. 

1) The spectrum of the wideband time-domain 
spread-spectrum signal s a (t) is subdivided 
into overlapping segments; typically the fre- 
quency overlap required is 10-15 percent. 

2) The power in each segment is normalized for 
equality. This lookup table of normalization 
constants is also known by the receiver to al- 
low renormalization in the reconstruction 
process. This normalization procedure ensures 
constant SNR for each received segment. The 
inverse Fourier transform of the normalized 
segment spectral data is the time-domain seg- 
ment data transmitted. 

3) Each time-domain segment is transmitted se- 
quentially with the appropriate RF center fre- 
quency of the Wi-Fi channel. Additionally, 
there is a time period between each transmis- 
sion to allow the receiver to retune and stabi- 
lize the local oscillator (in phase) for reception 
of the next transmission. 

4) The receiver calculates the spectrum of each 
time-domain signal segment, and uses the 
normalization table to correct the received 
spectrum of each segment to generate the 
wideband signal. 

5) The relative phase between adjacent segments 
is estimated using the procedure described in 
detail in Section 5.2 following. 

6) The differential phase data is used to rotate the 
phase of each segment to reconstruct the (co- 
herent) wideband spread-spectrum sig- 
nal S 'x(f) ' 

7) The despreading of the wideband signal is per- 
formed using the computationally efficient 
Fourier method. In particular the correlation 
diagram (correlogram) is generated by the op- 
eration 

C(r) = FFT^FFTlU^FFTMf))*) (6) 


Note that the Fourier transforms of the received seg- 
ment signals required in step (4) are used in the calcu- 
lation of the correlogram, so little additional computa- 
tions are required to calculate the differential phases. 
Also note that the Fourier transform of the transmitted 
signal is invariant, and thus the results are stored as a 
data table in the receiver which uses the wideband 
time-domain correlogram to estimate the TOA of the 
signal. 


Recovering the Differential Phase 

Consider two segments of overlapping complex spec- 
tral data. The receiver output can be represented as 

Segment 1 : 

Segment2 : A 2 S(f)e j ‘ />2 


where the (random) phases ^ are unknown, but 
the relative amplitudes are known by the receiver 
from the normalization table described previously. To 
recover the differential phase while reducing the ef- 
fects of noise, the following frequency-domain correla- 
tion is performed 

C - j (as (f)e* )(A 2 Sif)e j * ) df = Ae j ^E seg 
0 (8a) 

A/ 

where A = A 1 A 2/ E seg = j|s(/)|"d/is the energy in the 
o 

overlapping segment, and receiver noise has been ig- 
nored. As the segment energy is a real number, the 
differential phase can be recovered by 

<h “02= = ar g(c) I'fivn 


While the above analog correlation process can recover 
the differential phase, in practice the receiver output 
will include noise, and the signals are in a sampled 
form. Thus from (8a) with processing of sampled data 
the correlator output will be 


C-Ae 


M \-<h) ' 




/=i 


=Ae j ^~^'£ i C i = Ae 


ibh-if), 


C 


(=i 


(9a) 


where the ith sample of the digital noisy complex 
spectra is given by 

Segment 1: S[ = (s' R + n‘ n )+ y(s) +n‘ 12 ) ^ 

Segment 2: S 2 = (s« +n 21 )+ j(s) +n 22 ) 

where the random noise n has the statistical proper- 
ties anc j var [«] - (7 n ' ant | (-g e f our noise com- 


48 



Mobile Computing Volume 2 Issue 2, May 2013 


www.mc-journal.org 


ponents in (9b) are statistically independent. Unlike 
the integration in (8a), the summation in (9a) is com- 
plex, and thus performing (8b) using (9a) will result in 
an error in the differential phase estimate. Now con- 
sider properties of the real and imaginary components 
S ' 

of . The expected value of the real component of ith 

sample is £ M^I= ( 5 «) + ( 5 /) =P s (the signal power 
of the ith sample), so the summation over 

the N°vr samples in the overlapping region is propor- 
tional to the energy in the overlapping segment as in 
(8a). The corresponding variance is 


r Re(c,- 1 = E Re(c, 2 ) - E Re(c,- ) ~ 


= 2er„ 


+ G n 


( 10 ) 


Since the transmitted signals are pseudo-random, the 
real and imaginary spectral components are quasi- 
random and statistically independent, with means 


£[5 fi ] = £[5 / ]=0 


and variances 


r[.S' A .]= V ar[.S , / ] = cr s 2 


Applying these conditions to the 
samples in the segment gives 


N„ 


overlapping 


i[Re(c)]= 2N ovr of 

'ar[Re(c j = 4N ovr crj c 2 (l + I / 2y s ) 


( 11 ) 


where the pre-correlation segment spectrum SNR is 


2 / 2 
Ys = !a n 


, as the (normalized) segment power 


is A(7s and the noise power is 

(r) 


2 a 
put SNR 

r = 


2 at 


for the real component of is 


From (11) the out- 
C 


(2 N mr a;f 


<J~N 

s ovr 


N„ 


4N ovr cj;a 2 n {] + ]/2r s ) cr 2 (l + 1 / 2y s ) \ + \!2y s 


(12) 

Therefore this correlation has a process gain of 
■^ovr/(l + l/2y s ) ~ N ovr ' j ag j. approximation is 

a consequence of the segment pre-correlation SNR 
typically being much greater than unity. A similar 


analysis for the imaginary component of 

£|lm(c)]=0 


C 


var Im(c) = 4V ovr cr;cr 2 (l + 1/2 7 J « 4 N a 


gives 


? 2 

,0% <T„ 


(13) 


Q 

The imaginary component of is responsible for er- 
rors in estimating the differential phase as defined by 
(8b). As the phase error in radians is approximately 
the imaginary component divided by the real compo- 
nent (for small phase errors), combining the first equa- 


tion in (11) with the second equation in (13) gives the 
STD of the phase error (in radians) approximately as 

£t[A(Z1] ~ ; a n (l + 1 / 2 / ; ) _ 1 1 + 1 / 2/ ~ 1 

2 N ovr cr 2 s \ y s JnZ (14) 

Q 

In deriving (14), as the mean of real part of is much 
greater than its STD, the small variation in the real 
part is ignored. Thus given a particular segment SNR, 
the STD of the differential phase error decreases in- 
versely as the square root of the length of the overlap- 
ping segment. However, increasing the overlapped 
part of the segment reduces the transmission efficiency, 
so in practice this is limited to a small fraction of the 
total segment length. Consider a typical numerical 
example, with parameters similar to the WASP system. 
To ensure adequate performance the pre-correlation 

signal SNR is ,/s ~ ^ dB, and the length of the overlap 
is N°vr ~ 64 sam p] es ou t of a total segment length of 

512 samples. Applying (14) gives radians. 

While this phase error is quite small, the variance of 
the differential phase error (relative to the first seg- 
ment) increases proportional to the segment number, 
so the relative phase error in the last segment can be 
quite large. The cumulative effect of these errors is 
considered in the next section. 

Reconstructed Signal Correlation Performance 

Having determined the statistics of the differential 
phase error between adjoining segments, the overall 
effect on the wideband time-domain correlation of the 
spread-spectrum signal is now estimated. Although 
the wideband spectrum power varies across the band 
(maximum in the center and approaching zero at the 
edges), the power in each transmitted segment is nor- 
malized so that the transmitted power in each segment 
is similar, and hence the SNR in each segment is ap- 
proximately the same. As a consequence, from (14) the 
STD of the phase error in each segment will be similar. 
Further, the contribution of each segment to the over- 
all time-domain correlation (without phase errors) will 
also be approximately similar. Thus using the first 
segment as a phase reference, the relative effect of the 
differential phase errors on the correlator output is 
given by 

T e =1+ X ( 15a ) 

S=1 

where 
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® S =X A ^> ( 15b ) 

H=1 


is the summation of phase errors accumulated in a 
random-walk fashion across the reconstructed re- 
ceived time-domain signal. However, the magnitude 
of (15a) would involve a square-root operation, mak- 
ing statistical calculations difficult. Thus as an alterna- 
tive, the statistics of the magnitude-squared (power) is 

N 

analysed, so that with seg segments the relative 
power is 



f N seg- 

z 

( N seg _ 1 A 

r =r r * = 

cc c c 

1+X 

/O . 

e J 5 

1 + Z^‘ 


l *=1 

y 

V s=l J 


(16) 


seg —i seg — z, seg — i 

= N seg + 2 ^cos(® s )+ 2 2^ £cos(<t s -(I),) 


5=1 


5=2 t=s + 1 


To further simplify the statistical calculations associ- 
ated with (16), as the phase terms are generally rela- 
tively small (see numerical calculation in Section 4.2) 
the cosines in (16) are approximated by means of the 

truncated Taylor series cos x ~ ^ x * “ . Thus the ex- 
pectation (mean) of (16) can be approximated by 


N ,eg-1 r 1 N ,e S -2 N scg-l r 1 

E[r cc ]=N seg 2 - 2 E[cs> s 2 ]- 2 2 4 s *-®*) 2 ] (17) 

5=1 5=2 1=5+1 


Equation (19) gives the loss in the process gain due to 
the phase errors. For example, from the example in 

Section 4.2 with ° ^ ^ , applying (19) gives a triv- 

ial loss of 0.13 dB. Fig. 3 shows the loss in Process Gain 
as a function of the number of samples in the overlap 
and the input SNR. For minimal impact on the per- 
formance, the number of samples should be greater 
than 50 and the input SNR greater than -3 dB, which is 
realized by the WASP system. 


Process Gain Loss 
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FIG. 3 VARIATION IN THE PROCESS GAIN LOSS AS A FUNC- 
TION OF THE NUMBER OF SAMPLES IN THE OVERLAP AND 
THE INPUT SNR (y s ) . 


By applying (15b) and assuming the differential phase 
errors of each segment pair have the same variance 
and are statistically independent (due to the random 
phase associated with the frequency scanning process 
by the receiver) the expectations in (17) become 


Z £ h 2 J= Z £ 


r s 

Z A ^. 

U=i 


= Z Z £ K 


N \N - 1 1 

1 y seg v seg / 


(18) 


2 "seg - 1 


x X4v®,) 2 j= Z Z £ 

5=2 f=5+l 5=2 f=5+l 

N seg seg ~^fj^seg ~ ^)"|_ 2 


( s Y 


Zm. 

V»=i ) 


where ° is given by (14). Substituting these results 
into (17) and normalizing by the zero phase 

N 2 

error value seg gives the expected relative correlator 
output power as 



= 1 - 


kY 1 )" 


6 N 


seg 




N 


seg*-' A</» 
~6 


(19) 


Frequency Offset Estimation 

The analysis in Section 3 assumed perfect frequency 
synchronization throughout the mesh network. This 
section describes a frequency synchronization tech- 
nique, and estimates the accuracy of the synchroniza- 
tion process as well as its effect on TOA measurements, 
for both stationary (anchor nodes) and moving nodes. 

Accurate frequency synchronisation in the mesh net- 
work is required to ensure that clock drifts between 
transmitting and receiving a reply do not result in RTT 
measurement errors. For example, if the two TOA 
measurements associated with the RTT measurement 
are one second apart in time, then a differential drift of 
1 ppb in the clocks represents a RTT measurement er- 
ror of 1 ns. As typical low-cost oscillators have fre- 
quency errors of the order of 1-10 ppm, the system 
must provide a frequency synchronisation correction 
to reduce this initial error by at least a factor of 1000. 

Fig. 4 shows the measured variation in differential fre- 
quency between two oscillators over a long period of 
time in a temperature-controlled (air-conditioned) of- 
fice environment. It is assumed in the following analy- 
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sis that over short periods of time (up to a few minutes) 
the differential frequency is considered constant. 
However, over longer periods of time the variation is 
of the order of 20 ppb, so that differential frequency 
tracking is essential if an overall differential frequency 
error not greater than 1 ppb is to be met. 


Frequency Error Parameter 



FIG. 4 MEASURED DIFFERENTIAL FREQUENCY PARAMETER 
BETWEEN TWO OSCILLATORS OF NOMINAL ACCURACY OF 
1 PPM. THE FREQUENCY DIFFERENTIAL OVER THE 15 HOURS 
VARIES BY ABOUT ±10 PPB. THE DATA IN THIS GRAPH HAVE 
BEEN SMOOTHED TO REMOVE SHORT-TERM MEASUREMENT 
NOISE. 

Frequency Synchronization Procedure 

The basic concept to estimate the relative drift between 
two oscillators has been described in (Yu, Guo & Hed- 
ley 2009), but the following analysis provides more 
information on the process as well as analysis of the 
stationary case, and describes a new method for the 
practically important moving node case. 

The basic method to estimate the differential fre- 
quency is to determine the rate of change in the TOA 
between measurements separated in time. The TOA 
estimate from (3) can be modified to account for a lin- 
ear drift in the clock over a time period (t) as measured 
by the local clock. Defining the linear clock phase 
function as 

Off, a) = mod(^° + at,T pn ) (20a) 

the resulting TOA estimate as a function of time is 
given by 

TOA mn (f) = Off, a m ) + (/?„,„ + 5R mn )/c + A„ - Off, a n ) + T noise 
(20b) 

where at t=0 the clock phases are their nominal val- 

,o 

ues ^ based on their local clock, and the phases are 
assumed to change linearly over subsequent time. The 
time (t) in (20a) for the frequency offset estimation will 
be over many pn-code periods up to (say) one minute. 


i0 

and thus it is implied from (20a) that ^ is a constant at 
the start of transmissions over many pn-code periods 
if there is no frequency error; in practice for the 
method to work the transmissions must always com- 
mence at the start of a pn-code period, or 

mod (t,T pn ) 0, w ^ crc t measurec | by (-he local clock. 
Additionally in (20b), a range error has been added 
due to multipath propagation, plus a random meas- 
T 

urement noise no,se (with a zero mean). The mod 
function in (20a) complicates analysis, as the phase 
function has a sawtooth shape (Yu, Guo & Hedley 

T 

2009) with a sudden change of pn at a period 
T la 

of . As 01 is a small number, this sudden change 

in the TOA will occur infrequently, and as this change 
is of known value, the measured TOA data as a func- 
tion of time can be easily corrected by "unwrapping" 
the raw measurements (either up or down depending 
on the sign of the drift). With this correction, the mod 
operation in (20a) can be ignored in the analysis of 
performance, so if the two nodes are stationary and 
the multipath environment is static, then (20b) can be 
written in the form of a linear equation over time with 
additive noise, namely 

TOAnn ) = Pmn + a mn f + T noise ( 2 1 ) 

where a,nn ~ a,n ~ a " 

and & mn ~ ^ >m ~ ^ « ) + ( R " m + 8 R mn) /c + , a cons t a nt, at 

least over the frequency estimation period v' which is 
typically of the order of 10 seconds. If a series of TOA 
measurements is performed over time (typically with 
a constant update period of (say) 1 second), then a lin- 
ear regression analysis can be performed to determine 

the differential frequency drift parameter a " m . Note 

that the value of the constant parameter @ is not re- 
quired in determining the differential frequency pa- 
rameter, so that range measurement errors due to mul- 
tipath do not affect the estimation of differential fre- 
quency. However, the random component will result 
in errors in the estimate of the differential frequency 
parameter; and this uncertainty can be calculated us- 
ing the estimated STD in the TOA as given by (5) in 
Section 3.2. After unwrapping, the regression analysis 
on (21) results in an estimate of the STD in the differ- 
ential frequency parameter as 
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where <Te , the STD of the TOA measurement, is given 

T M 

by (5), is the total period of the sampling, and is 

the number of measurements ^ >> ^ . Applying (5) 
with a (low) SNR of 20 dB, and 1 1 measurements with 
an update period of 1 second (so T=10 seconds), (22) 
results in an (upper) estimated STD in the differential 
frequency parameter of 0.15 ns per second (or 0.15 
ppb). Thus if the nodes are stationary the frequency 
offset can be readily estimated to a precision that re- 
sults in trivial errors in the TOA, even if there are mul- 
tipath errors in the TOA measurements. 

The above procedure is satisfactory for static nodes 
where the changes in the TOA as a function of time 
can be interpreted as a consequence of a differential 
frequency between the two nodes. However, if one 
node is mobile, the changes in TOA can also be due to 
movement of the mobile node so that in this case the 
measurement of the differential frequency by the 
above method is not possible. 


Estimation of Range Error due to Motion 

Before analysing the problem of determining the fre- 
quency offset parameter for a moving node, consider 
the range estimation performance for moving nodes 
without frequency offsets. The RTT measurement pro- 
cedure consists of the mobile node transmitting and 
stationary nodes within radio range replying. Each 
node in the mesh network has an assigned time slot 
for transmissions, which incorporates both the spread- 
spectrum signals used for TOA measurement and data 
transmission. For example, the WASP system in (Sharp, 

Yu & Sathyan 2012) has ^ s ~ ^ time slots of 10 ms 

(t ) (t ) 

period v s > , so the position update period v «' is 0.6 

second. The delay in the reply depends on the time 

slots assigned to the mobile unit and the replying 

nodes, so a delay of p time slot periods will be in the 

T T 

range s " . From (4) and ignoring any frequency 
synchronization error, the RTT will be 


RTT m,n=~ 


( A ffl + A J- 


( 2 ^m« +AR »»i )o + (fl+dR,,,,, )„ 


\ + T„, 


(23) 

where D is the effective differential distance due to 


movement by the mobile over the interval from the 
time of transmitting to receiving a reply, and is given 
by D = \pT s cos 0 , where v is the speed, p is the slot 
delay parameter, and 0 is the angle between the veloc- 


ity and range vectors. For the WASP example above 
the maximum distance travelled in the update period 
is 0.9 m, where it is assumed that the mobile unit is 
attached to a person walking at a speed of 1.5 m/s. As 
the actual distance will vary dependent on the geome- 
try and slot numbers, this distance can be considered 
as quasi-random, so that from (23) the range error due 
to motion is given by 

■A = \ [ypT s cos 0 + (SR 0 + SR )]+ cT i 

2 (24) 

which has three statistically independent components, 
one associated with the movement of the mobile node, 
the second associated with multipath propagation, 
and the third associated with receiver noise. Now con- 
sider the statistical properties of the motion-related 
random process. The slot delay parameter p can be 
considered as having a statistical uniform distribution 
[ N s ! 2, N J 2] w j a mean d e l a y G f ^ ~ , and the an- 
gle between velocity and range vectors can be consid- 
ered as statistical uniform distribution over t n,n \ : 
where it is assumed that the replying nodes "sur- 
round" the mobile node in a statistically uniform 
manner. Thus averaging over all anchor nodes and 
positions in the coverage area the expected value of 
the range error associated with motion is 

E[s v ] = — e [p cos d] = — e[ p ]e[cos d] = 0 

2 2 (25 a ) 

as both expectations in the last expression in (25a) are 
zero. Therefore on average the motion of the mobile 
does not result in any range error measured at a point 

in time delayed by “ relative to the initial trans- 
mission by the mobile node. As the unbiased range 
error reference time will vary with each node pair, for 
maximum accuracy in position determination a veloc- 
ity vector correction should be applied. In most cases 
this complication is unnecessary. Similarly the vari- 
ance of the motion-related error can be calculated to be 

var k ] = (^) E [p 2 Wosef\=\ 

(25b) 

For example, when using (25b), the numerical values 
in the above WASP example give the STD in the range 
as 0.092 m. In an indoor environment the random mul- 
tipath component of the range error STD (for the 
WASP system data, see Table 1) is typically of the or- 
der of 1 m for static range measurements averaged 
over a large coverage area (Sharp, Yu & Hedley 2012, 


v Ja_ 

2 


N 

1 Y S 

12 


K) 

96 


52 



Mobile Computing Volume 2 Issue 2, May 2013 


www.mc-journal.org 


Sharp, Yu & Sathyan 2012). However, if the mobile 
unit moves a distance greater than about a wavelength 
in a NLOS propagation environment, the two TOA 
errors associated with the RTT measurement will be 
largely uncorrelated due to the random Rayleigh scat- 
tering of the radio signal. For the WASP case with a 
wavelength of about 5 cm, and for the above walking 
case, the distance travelled in one slot period is 1.5 cm, 
so nearly all of the slot delays will result in uncorre- 
lated measurements. Thus from (23) with two statisti- 

cally independent random variables ° , the effect of 
the averaging two TOA measurements reduces the 

-Ji 

range error STD by about a factor of " . Note that this 
applies only to the random component, but bias errors 
will be essentially constant over small distances, and 
thus the frequency synchronization process is not af- 
fected. However, it is shown in (Sharp, Yu & Sathyan 
2012) that mounting a mobile unit on the body of a 
person causes shadowing effects which results in in- 
creased GDOP, so the overall positional accuracy 
(product of GDOP and ranging error STD) may be 
worse than that with stationary base stations, despite 
the improved ranging error STD. 

Frequency Synchronization for Mobile Nodes 

While the frequency synchronization described in Sec- 
tion 5.1 is a simple procedure involving just TOA 
measurements between two stationary nodes, a modi- 
fied method is required when one of the nodes is mov- 
ing. One possible method to determine this motion 
status is by means of processing the data from an ac- 
celerometer in the mobile node. In particular, if the 
measured acceleration "noise" is below a threshold, 
then the node is deemed to be stationary, and the dif- 
ferential frequency parameter is calculated as de- 
scribed in Section 5.1. When the accelerometer "noise" 
is above the threshold, a previously measured (when 
stationary) differential frequency parameter value 
could be used. However, from Fig. 4, even under ideal 
conditions, the differential frequency parameter will 
vary slowly. If (say) the maximum allowable drift is 
1 ppb, it is suggested from Fig. 4 that the period in 
which the previously estimated value can be used is at 
most several minutes. The length of time will depend 
on the quality of the oscillator, so if this method is to 
be used, good quality (and expensive) oscillators are 
required. As the design goal is for cheap hardware, 
such a strategy is not ideal due to its unreliability 
when the period of motion without stopping exceeds a 
few minutes. Thus an alternative method of estimating 


the frequency offset parameter in nodes which are in 
motion is highly desirable. 

Before describing the method, it is useful to further 
refine the method of defining the offset frequency pa- 
rameter in each node. In the method described in Sec- 
tion 5.1 only a differential parameter between two 
nodes is required, but for the moving node case, the 
estimation of the frequency offset parameter for all 
neighbouring nodes in the mesh network is required, 
which will be explained below. As only relative fre- 
quency synchronization is required, not absolute syn- 
chronization, it is appropriate for one node in the 
mesh network to act as the network frequency refer- 
ence, and all other nodes are synchronized to this ref- 
erence node. By definition, this node now is deemed to 
have a zero offset frequency parameter, 
cc — 0 

namely " :/ . Further, as all stationary nodes which 

communicate with the reference node can determine 
the differential offset frequency relative to the refer- 
ence node, these nodes also can determine their fre- 
quency offset parameter, which can be transmitted in 
communications with other nodes. Thus by this pro- 
cedure, with communications throughout the mesh 
network, all stationary nodes (but as yet not moving 
nodes) can update in real-time their frequency offset 
parameter, and transmit this information to other 
nearby nodes. This procedure will result in a slow in- 
crease in the frequency offset errors as the number of 
hops from the reference node increases, so that 

N, 

with op hops, the STD in the frequency offset pa- 

, -|| , JN hop < - r ( c 0 , ■ ■ 1 

rameter will be v , where v is given by 

(22). However, this reduction in the accuracy can be 
readily compensated for by increasing the measure- 
ment time T in (22), so the accuracy estimate of 0.15 
ppb in Section 5.1 can be maintained in the multi-hop 
case. Furthermore, the frequency offset correction 
process used in correcting raw RTT measurements still 
only uses the relative offset between adjacent pairs of 
nodes, so the range error analysis in Section 5.2 re- 
mains valid for multi-hop stationary nodes. 

Having established the frequency offset parameter in 
each stationary node, the next requirement is to estab- 
lish this parameter for mobile nodes in real-time 
whilst in motion. The proposed method is similar to 
position determination using pseudorange data, 
whereby a range offset common to all measurements 

is determined as well as the ^ ) position. Because of 
this similarity with pseudorange least-squares (LS) 
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positioning described elsewhere (Yu, Sharp & Guo 
2009), only a summary of the details will be given here. 
From (4) and (23), the range measurement between a 

mobile ( m ' and stationary node &) can be expressed 
as 


D _ 

Ix mn 


cRTT „ , 


f k'D 1 
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« ( ? 2~ h ) 


<5R ni „+ SR . 
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mn nm 


d„ 


K+ A„ 


■ + cT„, 


ca, 


i ( f 2 h ) 


(26) 


where r ' is the time of the transmission from the mo- 
bile to the stationary node, ?2 is the time of the reply, 

and mn is the effective distance as defined in Section 
5.2. The first term in square brackets in (26) is either a 
measurement or a known value, while the second term 
is only a function of the mobile node (independent of 
the stationary nodes) with the only unknown the 
sought-after mobile node frequency offset parame- 
ter a ' n , and the last term is a random variable, which 
in Section 5.2 is shown to be dominated by the range 
errors associated with the indoor multipath environ- 
ment. 


The scattering of radio signals indoors is complex and 
depends on the details of the indoor environment, so 
theoretical calculations are difficult. However, experi- 
mental measurements (Alavi & Pahlavan 2006, Bellusci 
et al 2008, Gentile & Kik 2006, Alsindi et al 2009) and 
WASP range error data (Sharp, Yu & Hedley 2012, 
Sharp, Yu & Sathyan 2012) show that indoor NLOS 
ranging errors have a bias error which increases ap- 
proximately linearly with range, and a random com- 
ponent which is largely independent of range, pro- 
vided a small percentage (typically less than 10 per- 
cent) of large range errors are excluded; this exclusion 
is justified as they will be rejected in the position de- 
termination process (Yu, Sharp & Guo 2009, Sathyan et 
al 2011). Thus the indoors NLOS range error model 
assumed for the following theoretical analysis is 

AR = (SR 0 + XR)+ e r = P R + t' R ^ 7 ) 

where the term in brackets is the linear bias model 

and £r is a random component with zero mean. 
Thus if a LS analysis over all nodes in radio range of 
the mobile node is performed to determine the four 

unknowns the last three of which in this 

case are nuisance variables, then an estimate of the 


frequency offset parameter can be obtained. 


Using model (27) and (25b), (26), an analytical LS 
analysis (similar to that in (Yu, Sharp & Guo 2009, 
Sharp, Yu & Hedley 2012) shows that E[a m \ = 0 , and 
the STD in the frequency offset parameter of the mo- 
bile node is approximately given by 


A a m) f 


N r+ 5/9N r cT u 


1 24 

M 

In r 

l cT u) 


(28) 


_ 2 a; X^ mm )\{yT u f 

a AR ~ 1 77 > 1 777 AT 

where z 18 ^ , R is the number of 

stationary nodes in range, and max is the maximum 
radio range. In (28) it is assumed that the range meas- 
urements to the stationary nodes are statistically inde- 
pendent (because of the random indoor NLOS scatter- 
ing), and that these nodes "surround" the mobile node 
in a statistically uniform manner, as described in 
(Sharp, Yu & Hedley 2012). Notice that the STD is ap- 
proximately inversely proportional to the square-root 
of the number of nodes in range, and also inversely 
proportional to the update period; these relationships 
are similar to the stationary case given by (22). 

The simulation results of the frequency offset estima- 
tion process are shown in Fig. 5. The simulation is 
based on a mobile unit moving along a circular path 
with parameters as defined in the caption text for 
Fig. 5. 


Estimated Frequency Offset Parameter 



FIG. 5 SIMULATED PERFORMANCE OF ESTIMATING THE FRE- 
QUENCY OFFSET PARAMETER FOR A NODE MOVING AT A 
SPEED OF 1.5 M/S ON A CIRCULAR PATH OF RADIUS 25 M 
WITH AN UPDATE PERIOD OF 0.6 SECONDS. THE COVERAGE 
AREA IS 100 M SQUARE, WITH 44 NODES UNIFORMLY RAN- 
DOMLY DISTRIBUTED (BUT WITH A MINIMUM SEPARATION 
OF 7.5 M), OF WHICH (ON AVERAGE) 18 ARE WITHIN RADIO 
RANGE. THE RANGE ERROR MODEL PARAMETERS 
ARE A, = 0.03, a= 1 m, WITH A MAXIMUM RADIO RANGE 
OF 40 M. 
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The frequency offset is represented by a mean of 
5 ppm plus a sinusoidal component with amplitude 
200 ppb and a period of 200 seconds, similar to those 
assumed in the simulation in (Yu, Guo & Hedley 2009). 
This variation is not intended to model actual fre- 
quency variations, but to demonstrate the ability to 
dynamically track quite rapid variations. The mesh 
network parameters chosen are typical of the deploy- 
ment of the WASP system in an indoor office-type en- 
vironment. Fig. 5 shows both the raw estimate of the 
frequency offset parameter, as well as the Kalman fil- 
tered smoothed version. The raw data has a mean er- 
ror in the frequency offset of 1.6 ppb with a STD of 
8.1 ppb, compared with the analytical estimate from 
(28) of 7.8 ppb. The filtering is based on a Singer g-h-k 
Kalman filter (Brookner 1998) for tracking a target 
whose dynamics have a random acceleration with a 
exponential autocorrelation function with a time con- 
stant of 200 seconds. The STD of the Kalman filtered 
data is 4.8 ppb. As described in Section 5.2, the RTT 
reply period (T) is a pseudo-random variable (uniform 
distribution in the range 0 to 0.6 second in this case), 
so the associated STD in the range 
. ca(a V / Jl2 = 0.26 . 

is v ' v m. Assuming that the range bias 

errors are corrected for, so that only the random range 

errors ( a '' ~ ^ m in this case) remain, then the random 
effects in the estimation of the frequency offset pa- 
rameter will result in the effective STD in the random 

• , Vl 2 +0.26 2 =1.03 ryr, ., 

range error increasing to m. thus the 

effect on ranging accuracy by errors associated with 

the determination of the frequency offset parameter 

for a moving node is insignificant 

Measured Performance Results 

The methods described in the previous sections have 
been implemented in the WASP positioning system, so 
that testing of the theory is possible through meas- 
urements. However, as the overall performance is the 
integrated sum of the individual subsystems, such as 
the signal concatenation and frequency synchroniza- 
tion processes which are the focus of the methods and 
analysis described in this paper, the performance of 
these individual subsystems cannot be analysed in 
details. However, it is confirmed that the ranging per- 
formance meeting expectations is indicative that the 
subsystems also perform within specifications. 

The WASP performance data are summarized in Ta- 


ble 1, and are listed in chronological order. The data 
include both range (Ar) and position errors (e), mean 
and median statistics, static and dynamic measure- 
ments, and LOS/NLOS propagation environments. 
The data for the first four cases in Table 1 are in the 
indicated referenced paper, while the final case is re- 
ported here for the first time. Note also that the num- 
ber of nodes increases from the initial 9 to the final 46 
nodes; and that the corresponding area covered in- 
creases from 225 m2 to 10,000 m2. The final example is 
for the WASP system deployed in a complex of build- 
ings, much as would be expected in a commercial de- 
ployment. Although there are a wide variety of data 
presented in Table 1, the general conclusion is that 
LOS performance (range and position) is 15 cm or bet- 
ter, while NLOS performance (in a building) is gener- 
ally better than 1 m. 

The static LOS measurement errors at high SNR (about 
50 dB) are dominated by errors associated with the 
wideband reconstruction and frequency synchroniza- 
tion subsystems. The 15 cm STD (or 0.5 ns) range 
measurement noise can be attributed to "reconstruc- 
tion noise" associated with the wideband signal con- 
catenation process, and as such represents the lower 
limit in possible ranging performance. This lower limit 
in ranging also places a lower limit on the accuracy of 
the frequency synchronization process. Applying (22) 
with M=ll and T=10 seconds yields a trivial frequency 
error STD of 0.055 ppb. Even at the typical lower oper- 
ating SNR of 30 dB, the estimate frequency STD would 
be 0.55 ppb. In these cases, it is clear that frequency 
synchronization errors will not be a problem. The 
NLOS measurements listed in Table 1 have ranging 
errors of the order of 0.5-1 m. In this case, the indoor 
multipath propagation environment is the dominant 
factor, rather than the performance of the analog radio 
and its associated signal processing. These WASP per- 
formance data are much better than other reported 
results (as summarized in Section 2) by at least a factor 
of 5, which is about the factor (125 MHz to the 20 MHz 
Wi-fi channel ratio) of the wideband reconstruction 
process. This result suggests that the reconstruction 
process described in Section 4 is effective in a multi- 
path environment, with the improvement factor close 
to the bandwidth enhancement factor. Ultimately, per- 
formance in the NLOS indoor environment is limited 
by available bandwidth in the ISM bands, rather than 
the signal reconstruction process. 
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TABLE 1. SUMMARY OF PERFORMANCE DATA OF WASP SYSTEM AS REPORTED IN THE INDICATED PAPER. 


Paper 

LOS 

^ r rms 

NLOS 

A r 
rms 

Nodes 

Area, d 

Static NLOS 

P 

° pos 

Static NLOS 

£ pos 

Dyn NLOS 

£ pos 

Dyn LOS 

£ pos 

Hedley et al, 2008 

15 cm 


9 

225 m 2 , 5.0 m 


49 cm 



Sathyan et al, 2011 

15 cm 


9 

225 m 2 , 5 m 

48-61 cm 



7 cm (3) 

Sharp, Yu & Sathyan, 2012 



29 

2100 m 2 , 8.5 m 

48 cm 


76 cm 


Sharp, Yu & Hedley, 2012 


108 cm 

36 

2100 m 2 , 7.6 m 

46 cm (1),70-100 cm (2) 






94 cm 

46 

10000 m 2 , 14.7 m 


65 cm 




(1) Core locations in building, (2) Positions at the edge of buildings, (3) Based on two mobile units 0.7 m apart on a bike in a velodrome - meas- 
ures differential position error only. The symbol signifies the median. 


Conclusion 

This paper investigated the application of IEEE 
802.11/Wi-Fi analog chip-radios to positioning systems. 
Because of the desire for cheap implementation in 
commercial applications, it is proposed that the Wi-Fi 
chip-radios used for data communications are also 
used for radiolocation applications in indoor envi- 
ronments. Achieving the desired sub-meter ranging 
accuracy with this type of radio is challenging, as the 
design of these chip radios is aimed at data communi- 
cations rather than radiolocation. To achieve the goal 
of accurate radiolocation, two main areas where inno- 
vative design is required have been identified, namely 
achieving an effective wide bandwidth, and frequency 
synchronization in a mesh network with mobile nodes, 
all in a radio-hostile indoor multipath propagation 
environment. By using commercially available analog 
chip-radios and appropriate signal processing tech- 
niques, it has been shown that indoor ranging accu- 
racy of about 1 m can be obtained, and that when ap- 
plied in a mesh network, positioning accuracy of 0.5 m 
can be obtained without the need to adopt ultra-wide 
bandwidths and the associated high frequency (and 
power hungry) electronics, or complicated timing syn- 
chronization schemes. Current generation Wi-Fi chips 
have close integration of both the analog radio and 
digital signal processing functions into a single chip, 
so that the methods used in the WASP system could 
not be readily adopted with such chips. However, it is 
suggested that by adding further specialized signal 
processing hardware, a chip with both Wi-Fi data 
communications and accurate radio location function- 
ality could be developed without much difficulty. 
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